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Abstract Dense Bose superfluids, as Hell, differ from dilute ones by the 
existence of a roton minimum in their excitation spectrum. It is known that 
this roton minimum is qualitatively responsible for density oscillations close 
to any singularity, such as vortex cores, or close to solid boundaries. We show 
that the period of these oscillations, and their exponential decrease with the 
distance to the singularity, are fully determined by the position and the width 
of the roton minimum. Only an overall amplitude factor and a phase shift are 
shown to depend on the details of the interaction potential. Reciprocally, it 
allows for determining the characteristics of this roton minimum from static 
"observations" of a disturbed ground state, in cases where the dynamics is not 
easily accessible. We focus on the vortex example. Our analysis further shows 
why the energy of these oscillations is negligible compared to the kinetic 
energy, which limits their influence on the vortex dynamics, except for high 
curvatures. 

Keywords Quantum fluids and solids • Superfluid phase • Hydrodynamic 
aspects of superfluidity 

PACS 67. • 67.25.D- ■ 47.37.+q 



1 Introduction 

Hell is the low temperature, low pressure superfluid phase of 4 He. The ro- 
ton minimum, in its excitation spectrum, has been inferred by Landau [1] 
from the viscosity measurements of Andronikashvili [2]. It has been shown, 
by Feynman [3] to be due to the dense packing of 4 He atoms. Solidifica- 
tion, which occurs above 2.5MPa for low temperature 4 He, can be seen as a 
condensation of rotons, due to their interactions. 
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The low value of the excitation energy at the roton minimum also suggests 
that the superfiuid has a strong susceptibility for spatial perturbations of 
wave number k a , the position of the roton minimum. Localized perturbations 
will then produce oscillations in the superfiuid density in their neighborhood 
[5-7]. This has been known for long. For instance, using an ansatz for the 
superfiuid order parameter, based on the only neighborhood of the roton 
minimum, T. Regge [4] could show the existence of these oscillations close 
to the free surface of superfiuid 4 He. He could obtain the surface tension 
within 30% of the experimental value. More recent theoretical works [5,7] 
obtain such oscillations in the order parameter close to other singularities, 
as vortex cores. However, the sensitivity of these oscillations to details of the 
real problem at hand is not clear. 

Up to now, because of the smallness of the scales involved in this phe- 
nomenon (e.g. of the order of the Angstrom in 4 He), any experimental inves- 
tigations remain a great challenge. For instance, one could think of neutron 
diffraction by a vortex lattice, created in a 4 He film by a rapid rotation. 
However, rotation frequencies as large as 10 5 rd/s would result in a distance 
between vortices of order 1 /im, which ask for resolving orders up to 10 3 in 
the reciprocal lattice 

So then, two approaches can be followed. First "numerical experiments" 
through ab initio calculations [8] allow to compare the ground state and the 
one vortex state for a limited number of atoms. Secondly, Gross-Pitaevskii 
equations are often proposed for modeling the dynamics, and space depen- 
dence of the order parameter. The interaction term in these equations is 
chosen so as to fit the dispersion curve of elementary excitations. But a large 
class of such interaction terms can fit the same dispersion curve. What is the 
influence of their differences? What is the influence of details of the disper- 
sion curve such as the phonon (long wavelength), maxons (maximum of the 
dispersion curve) or very short wavelength regions? Up to what precision do 
we have to fit the roton minimum (simple parabola, skewness)? At the end, it 
is known that, for strongly interacting superfluids as 4 He, density could not 
be simply proportional to the squared order parameter, as in dilute systems. 

Determining the parameters which control the extension and amplitude of 
these oscillations, and discussing their implications are the goals of this paper. 
We shall focus on density oscillations far from the singularity, where they are 
small and amenable to a linear equation. We shall compare our approximation 
both with "exact" resolution of various Gross-Pitaevskii models, and with the 
results of ab initio calculations. 

The existence of density oscillations close to a vortex core is expected to 
have tremendous importance in the study of quantum turbulence in which a 
central role is played by vortex reconnections [9,10] and thus it is useful to 
quantify how these oscillations depend on the shape of the roton minimum. 
In particular, this study could be necessary to the design of numerical inves- 
tigations of vortex reconnections, as it was done formerly without taking into 
account the roton gap [11]. It is why we shall focus on the vortex problem. 

This paper is organized as follows. First and foremost, in section 2, the 
universal influence of the roton minimum close to a localized perturbation 
is shown on a simple example, the Pomeau-Rica version [6] of the Gross- 
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Pitaevskii equation [12]. In the following section 3, we extend this linear 
approximation to density oscillations in real superfluids. The order parameter 
oscillations around a vortex core are then compared, section 4, in the Pomeau- 
Rica model, to our linear approximation. Then, in section 5, we apply the 
same comparison to a more elaborated Gross-Pitaevskii model [7], aiming 
at better reproducing the full dispersion curve than in the Pomcau-Rica 
version. Our last comparison, section 6, concerns the density profile around 
a vortex obtained via ab initio variational calculations [8] . The last section, 7, 
is devoted to energy considerations. We finally conclude on the implications 
for real systems. 



2 Linear treatment of the non-local Gross-Pitaevskii equation 

The Pomeau-Rica [6] version of the Gross-Pitaevskii [12] equation reads: 

ift^(r,t) = -|^(r,t)-/i^^ (1) 

where U o 0(\r — r'|/a) represents the interaction potential between two atoms 
located at r and r', a is the range of the potential. The interaction 0(x) is 1 
if |x| < 1, for |x| > 1 and /j, is the chemical potential, fixing the equilibrium 
density n = \ip(r,t)\ 2 . 

Using the transformation r — > ar, ip — > y/nip, t — > (2ma 2 /~h)t, Eq. 1 can 
be written in the following non-dimensional form: 

idt1>(T, t) = -A^(r, t)-47r^(r, i)/3-M^(r, t) J d 3 r' 6(\r - r'|)^(r', t)\ 2 (2) 

where 

2ma 2 U na 3 

A = f? {i) 
remains the only parameter. Any static perturbation will appear as an ad- 
ditional term V(r) = -0(r)V(r) in the right hand side of equation (2). Far 
from the perturbation, the solution of Eq. 2 can be written to linear order 
as V = 1 + Sip, with dtp and its complex conjugate dip* verifying the linear 
system 

id t 5tp(r,t) + A5tp(r,t) 

-A J d 3 T'9(\r-T'\)(SiP(r',t) + Sr(r',t)) = V(T) (4) 
-id t 5tp*(r,t) + A5ip*(r,t) 

-A J d 3 r'8(\r-r'\)(5yj(r',t) + 6r(r\t)) = V(r) (5) 

By Fourier transforming, we get: 

{uj-k 2 - A6{k))8~ip{lt,Lj) - A0(k)Sip* (k,w) = V'(k) 
(w + k 2 + A6(k))S~tp*(k, u) + A9(k)S~iP(k, u) = -V (k) (6) 
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where 6{k) is the Fourier transform of the normalized interaction potential 
9: 

6{k) = ^(sin(fc) - fccos(fc)) (7) 
Eliminating S~ip from the former system (Eq. 6) gives: 

(lu 2 -uj(k) 2 )S~^(k,uj) = (uj + k 2 )V(k) (8) 
where w(fc) is the dispersion relation of the linear excitations, namely 



LU 



(k) 2 = k 2 (k 2 + 2A9(k)) . (9) 



The frequency to = oo(k) allows a non zero Sip(k, go) without any perturbation 
(V^(k) = 0). The static (uj = 0) solution of equation (8) is: 

('») 

In this model, Sip (k) represents half the Fourier transform of the relative 
density deformation: 

Finally, the form of the deformation Sp(r)/p Q in the physical space is 
given by the inverse Fourier transform 



^(k) = -^- (ID 



W = 7^ / -(k)e Jk - r rfk . (12) 



Po [ ' (2tt)3 J Po 

In the case of interest for us, uj{k) 2 has a deep minimum (the roton 
minimum) at k — k a that will clearly dominate the right hand side (RHS) 
of Eq. 10. Let us write in its neighborhood w(fc) 2 = fl 2 + c 2 (k — k ) 2 , where 
fifl is the energy roton gap and to* = TiQ/c 2 the "roton mass". 



3 Extension to real systems. Linear vortex 

Equation (11) is central in our paper. Let us present here our arguments 
for extending it to any real superfluid, as 4 He. We are interested in the 
deformation 5p(r) far from the perturbation. As Sp/p a is small, it should 
obey a linear equation: 

/A = V(r,t) (13) 

Po 

Due to space and time translational invariance, the Fourier transform of 
this equation is local: 

#(w,k)(w 2 -w(k) 2 )^(w,k) = VYw.k) (14) 

Po 
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The shape we gave to this equation simply expresses that deformations 
can propagate in the superfluid, in the absence of any perturbation, if u) = 
±u>(k). Our point is that F is smooth in the neighborhood of k = k , as well 
as V, if the perturbation is localized. For a static perturbation, we can then 
write, in the general case: 



For a plane localized perturbation, as a free surface, or a wall, we can 
write V(k)/F(k) = v(k x )S(k y )5(k z ), where 8 is the Dirac function, and the 
plane is perpendicular to the x axis. Because the roton minimum is assumed 
deep, i.e. fl 2 is assumed small, the RHS of Eq. 15 has a sharp maximum at 
k = k . Thus, the very precise shape of the perturbation is not important as 
long as it is localized, i.e. the Fourier transform v(k x ) is smooth, and we can 
take it as a constant v(k a ). The result: 

> - - V M • <w) 

is similar to the ansatz proposed by Regge [4]. 
Let us now turn to the linear vortex case. 

To model the perturbation induced by a vortex line, we consider a lo- 
calized perturbation that depends only on the distance r from the z-axis 1 
(cylindrical symmetry). In this case, we can write V(k)/F(k) = v(k±)8(k z ), 
where 8 is the Dirac function. Again, due to the sharpness of l/uj(k) 2 close 
to the roton minimum, the very precise shape of the perturbation is not im- 
portant as long as it is localized, i.e. the Fourier transform v(kj_) is smooth. 
In this spirit, we can Taylor expand the perturbation as v(k±) sa (v(k a ) + 

v'(k )(k± — k Q ))5(k z ). Using Eq. 16, we can approximate y-(k) as: 
<5p v(k ) + v'(k )(k ± -k ) 

Then, using the approximative form of |^(k) (Eq. 17) and Eq. 12 , we 
get Sp(r) = Sp(r) with 



S P( . 1 f°° d(k ) + v'(k )(k ± -k ) 



where Jo{x) = (l/2n) exp(za; cos </))d(j) is the Bessel function of the first 
kind. The shape of the perturbation (Eq. 18) can be accurately numerically 
evaluated and is found oscillatory (data not shown) . Given the approximative 
form used in Eq. 17, it is expected to be valid only when it is small, i.e. for 
r relatively large. In the integral, we can thus consider k±r as large, for the 



1 Doing so, we neglect the influence of the phase, which changes by 2iv when 
turning around the vortex. 
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interesting k± values. This allows to use the following asymptotic form of the 
Bessel function 

, , 2cos(x — 7r/4) , . 

J (x) » V 1 > . (19) 

V 27TX 

Using the asymptotic form of the Bessel function (Eq. 19), the k± integration 
entering in Eq. 18 can be performed in the complex plane, which gives: 

dp. . exp(— &ir) / , , k\ , \ , . 

-£(r) ~ a F \_ 1 ' cos k a r - tt/4 + + O (20) 

with ki = Q/c and: 

7 K 

(7 = 



(2tt)V2 c 2fc 1 
7 2 =f;(fc ) 2 + fc?5'(fc ) 2 

cos(0 o ) = ^ (21) 

7 

the sign of 7 being chosen such that — | < O < -|. 

Thanks to Eq. 20, we can see that the density profile close to the vortex 
behaves as an oscillatory function of wavelength 27r/fc D that exponentially 
tends to zero with a characteristic length scale 1/fci. Let us emphasize that 
the obtained shape of the deformation 8p(r) only depends of k r, k\jk and 
an additional phase-shift <j> . In the small ft limit (i.e. deep minimum), the 
ratio k\/k is expected to be small compared to unity. The sign and ampli- 
tude of 4> a depends on the shape and variation of the perturbation v(k±). 
If v(k a ) ~ —k v'(k ) as obtained with a Gaussian function of characteristic 
spatial extension the atomic size ~ l/k , we get a small negative phase-shift. 
If on the contrary, v(k a ) is small compared to k v'(k ), then (f> can take a 
significant value, which sign depends on the shape of v around k Q . We see 
that the qualitative behavior of the density profile close to a singularity is 
expected to be universal, i.e. independent on the precise shape of v, only an 
additional phase-shift is model dependent. 



4 Numerical solution of the Pomeau-Rica equation 

To compare the linear prediction obtained in Eq. 20 with the full rectilinear 
vortex solution of the Gross-Pitaevskii equation (Eq. 2), we perform a nu- 
merical simulation aimed at evaluating the ground-state of Eq. 2 under the 
constraint that the phase turns of 2tt along a line, as it was done in Rcfs. [12, 
5,7]. In this axisymmetric geometry, we are thus looking for an order param- 
eter of the form ip(r, t) = f(r)e lv , where (r, ip, z) is the cylindrical coordinate 
system. The amplitude f(r) is then given as the solution of the following 
boundary values problem 

r(r) + ^+(f™l-l)/(r)= (22) 
Af(r) / / / eUr 2 + r' 2 - 2rr> cos ip + z 2 ) r' f 2 (r')dr'dtpdz 

Jr'=0 J<p=Q Jz = -oc ^ ' 
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Fig. 1 Comparison between the "equilibrium" density around a vortex (dots) in 
the Pomeau-Rica model, with the first order approximation Eq. 23 (full line). Eq. 
22 is solved via a relaxation method. ki/k a = 0.0659, 4> = 0.1694. 



subject to the conditions /(0) = and /(+oo) = 1. 

Using the non-local potential 6 of Pomeau and Rica [6] , we can perform 
analytically the integration over z and (p using the Elliptic integrals of the 

first F{k) = Jo V1 _ fc 4yi-t* and second £ (*) = So dt kind ' We can 

show then that the RHS of Eq. 22 is given by, for r < 1 



8f{r)A lZ Vl-^-^f U y^^ J r'f(r')dr' 
+8f(r)A / -ff? [28 (b) - (a + r'f 2 {r')dr' 

JX-r 

and for r > 1 

/■r+l 

8/(r)yl / v^Vpf (6) - (a + 1)J"(&)] r'f{r')dr' , 

where a = r ^ r / -1 , & = — a. To numerically solve Eq. 22, we use a 

relaxation method [13]. The Laplacian is furthermore estimated using a joint 
Crank-Nicolson and Gauss-Seidel algorithms [13] to ensure numerical stabil- 
ity. Starting with the uniform amplitude f(r) — 1, the relaxation method 
converges towards the "equilibrium" distribution f(r) that solves Eq. 22. 

In figure (1), we compare this approximation with the "exact" solution 
of the Pomeau-Rica model, for A = 40, obtained with the above relaxation 
method. That is, we compare f(r) 2 with: 

p hn = 1 + —== cos(fc G r - - - 4>) exp(-fc!r) (23) 
V k r 4 
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a and <j) being fitted for the best agreement. The parameter A — 40, for 
which k a w 5.40, is close to the "spinodal" value A s ~ 43.43, for which 
J? = 0. As it can be seen in Fig. 1, the first order approximation Equation 23 
(solid line) almost perfectly fits the "equilibrium" density (dots). The value 
of ki/k is here k\/k D — 0.0659. The value obtained for <f) is <f> — 0.1694, and 
a = -0.9917. 



5 Melting pressure Helium 

In order to be closer to real superfluids, and to estimate to what extent the 
above ideas can apply to experimental systems, we consider now a more elab- 
orate model. It has been used in a series of papers by Berloff and Roberts [7] , 
to which we refer for details. Apart from a larger number of parameters in tai- 
loring the two-body interaction potential, it differs from the Gross-Pitaevskii 
approach by taking into account three-body interactions. Altogether, it al- 
lows to adapt the dispersion relation, both qualitatively and quantitatively, 
to that of the real superfluid of interest. 

In this model, the Schrodinger equation 1 is changed to: 

h 2 

ihd t ip(r, t) = -—AiP(r, t) - ^(r, t) 

+V(r,t) ^ cPv'U(v-v')\^{v',t)\ 2 + W\^(r,t)\ 4 ^j (24) 

where the three-body potential is treated as local, with intensity W. 

Defining a = h/yJ2mfi, the same transformation than above, r — > ar, 
V> — > \/nip, t — > (2ma 2 /h)t, yields: 

-idt1>(T, t) = Z\^(r, t) + ^(r, t) (1 - X |V>(r, t)\ 4 ) 

-^(r,t) i^j dVf/(|r-r'|)|^(r',t)| 2 ) (25) 

The reduced potential U is taken as (for convenience, we keep the same 
notations as [7]): 

U{r) = (a + fiA 2 r 2 + M 4 r 4 )e- AV + ^-^^ . (26) 

The dispersion relation is then given by: 

uj 2 = k 4 + 4fc 2 x + 2k 2 U(k) (27) 

where U(k) is the Fourier transform of U(r). 

The choice A = 1.3282, B = 0.1992, a = 11.5881, (i = -28.48, 8 = 
7.723, rj = 0.003775, x = 3.5, gives good agreement with the experimental 
dispersion curve at low temperature, at a pressure close to the melting one 
[14] (see figure 3). In particular, we obtain good values not only of the roton 
minimum and its curvature, but also of the velocity of sound and the "maxon" 
energy. 
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Fig. 2 Comparison between the "equilibrium" density around a vortex (dots) in 
the Roberts model Eq. 25, with the first order approximation Eq. 20 (full line). Eq. 
25 is solved via a relaxation method. ki/k ~ 0.1144, <j> = —0.99. 

Similar numerical investigations of the density profile near a vortex as 
the ones done for the Pomeau-Rica case (Eq. 22) can be performed. In this 
case, the RHS of Eq. 22 can be analytically computed using Bessel functions 
(see Ref. [7] for details) . The comparison between the "equilibrium" density 
and our first order approximation Eq. 20 is shown in figure 2. Again, the 
agreement is very good, the discrepancy being visible only close to the second 
minimum. The value of k\/k is here 0.1144. However, the value of <j) is 
4> = —0.99. Correspondingly, the hollow core of the vortex, that is the radius 
on which the density is nearly zero, is much wider with this model than with 
the Pomeau-Rica one. This quantitative phase-shift can be explained using 
a perturbation potential v that is close to around k Q . 



6 Comparison to variational calculations 

Up to now, we compared our linear approximation to Gross-Pitaevskii ap- 
proaches. The structure of a vortex has also been studied via variational 
calculation approaches using appropriate trial wave functions [8]. These cal- 
culations can be considered as numerical experiments on well controled fluids. 
In such a case, the authors have difficult access to the dynamical properties 
of their fluid model [15]. It is then interesting to test both if the shape of 
the density structure is compatible with our predictions, and if the so ex- 
tracted characteristics of the roton minimum match well the experimental 
measurements on the simulated state of the fluid. 

In figure 4, we compare the results of reference [8], corresponding to 4 He 
close to the melting pressure, to our formula for density variations close to 
a vortex line. The parameters are k a = 2.0710 m , and k\jk = 0.12, 
slightly different from those of the Roberts model treated above. However, 
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Fig. 3 Filled black circles: fit of the experimental 4 He dispersion curve (neutron 
scattering [14]). Filled red diamonds, and curve: parabolic fit, ki/k = 0.12, k = 
2.0710" 10 m _1 . Dashed line: Roberts' type model. 



the fit of the experimental [14] roton minimum with the formula: 

uj(k) 2 = c 2 {kl + {k-k f) (28) 

where c is here adjusted, is at least as good as with the Roberts' formula 
(see figure 3, dashed line). Several informations can be drawn from the above 
results. 

— First, the excitation energies corresponding to the trial wave function 
used in reference [8] present a roton minimum whose wavelength and 
mass closely correspond to the experimental ones. This is a very good 
test of the pertinence of this wavefunction, which is difficult to obtain by 
other ways [15]. 

— The phase shift is small = 0.0186. It is much smaller than obtained 
from a Roberts' model, even with a small three body term ((f> = 0.6650). 

— The amplitude of the density oscillations is smaller than with the corre- 
sponding Roberts' model. 

Let us comment on the last two remarks. The wave function of reference 
[8] has the peculiarity to give a finite value for the density on the vortex 
line, thus distributing the vorticity on a finite radius. No Gross-Pitaevskii 
model can do that. On the opposite, the Roberts type model tends to en- 
large the hollow zone around the vortex line, especially the three body term. 
It naturally shifts the oscillations outside, increasing (/>. It probably also en- 
hances the perturbation induced by the vortex, and thus the amplitude of 
the oscillations. 

This certainly raises a doubt about the utility of Gross-Pitaevskii ap- 
proaches. As we have seen, the general shape of the oscillations, i.e. their 
pseudo-period and their amplitude decrease, is completely determined by the 
characteristics of the roton minimum. All the remaining physics reflects in 
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Fig. 4 Comparison between the reference [8] density profile (black dots, from 
their figure 3, dotted line), and our formula 23 with k\/k = 0.12, a = —0.4643, 
= 0.0186. 

the amplitude and the phase of these oscillations. Unfortunately, the Gross- 
Pitaevskii models give model dependent results, even if they reproduce the 
whole dispersion curve. We thus cannot trust the amplitude and phase pre- 
dicted with these models. 

However, it is interesting to look at the dependence of the amplitude 
versus ki/k a , i.e. the non dimensional roton gap. In reference [8] the density 
profile around a vortex line has also been studied for low pressure 4 He, where 
ki/k ~ 0.16. The obtained amplitude of the first peak is very close to the 
high pressure case (not shown). This almost insensitivity of the amplitude 
to the roton gap appears also with the simplest non local Gross-Pitaevskii 
model, the Pomeau-Rica one, except may be very close to the zero gap, where 
the amplitude seems to fall rapidly (however, the time necessary to reach the 
equilibrium shape diverges in this limit). This agreement between the two 
very different approaches gives us confidence in the near insensitivity of the 
amplitude with the roton gap, which will be important in the next section. 



7 Energy considerations 

Even within our linear approximation of the density deformations, it is pos- 
sible to have a look at their contribution to the vortex energy. It comes from 
the fact that, in a free propagating linear excitation, the kinetic energy and 
the potential deformation energy are equal in average (this deformation en- 
ergy not only comes from the two body interaction potential between atoms, 
but also from the real part of the wavefunction gradient, the imaginary part 
giving the kinetic energy). Furthermore, mass conservation allows to relate 
the deformation amplitude of this excitation to the longitudinal velocity (i.e. 
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u(k) the component along k of the Fourier transform of velocity): 

dtp + djpuj = —u>(k) — (k) + ku(k) = (29) 

Po 



Thus the contribution of a wave vector k to the deformation energy is: 

2 

(30) 



2 fc 2 



5 Ak) 

Po 



and the total contribution of the oscillations to the energy per unit length of 
the vortex is then: 

6 = Jffis J e(k) d " k ■ (31) 

We are interested here in quantifying the energy induced by the oscilla- 
tions as given by the simplified form of the perturbation Eq. 17. In particular, 
we are not taking into account the energy deformation induced by the vortex 
core, this can be done by other means [12]. In this spirit, we restrict the do- 
main of integration in Eq. 31 to a neighborhood of the roton gap ko- As we 
saw, this approximation is consistent when the roton gap is deep, i.e. when 
ki = Q/c <C ko (see Eqs. 20 and 21). For example, the integration domain 
can be restricted to ko ± nki, with n an integer such that nk\ <C fco, and we 

get, using the approximated form of y-(k) (Eq. 17) 



2£i 

1 Po f k ° +nk > dk ± [v(k )+v'(k )(k ± -k f 



4?r fl 2 J k =feo _„ fel k± 



1 I ( k ± -k \ ' 



(32) 



In the integrand of the former integral (Eq. 32), the numerator can be ex- 
panded and gives three terms, two of them being proportional to (k± — k a ) 
and — fc ) 2 . Their contribution to the energy (Eq. 32) are of the order 
of n(fci/fco) 3 and, henceforth, will be neglected. It also shows that, in the 
limit of vanishing roton gap fci/fco — > 0, the phase-shift 4> contributes with 
a negligible factor to the energy. The remaining leading term, proportional 
to v(k ) 2 i=a 7 2 , 7 2 being itself proportional to a 1 (see Eq. 21), can be shown 
to be of the order of (fci/fco) and independent on n. We finally get for the 
deformation energy 

"»—2%—- ^ 

Looking at the expression of e (Eq. 33), we see that it contains k\ at the 
numerator. Let us remark that we checked numerically on the Pomcau-Rica 
model that the amplitude a 2 does not diverge when the parameter A tends 
to the spinodal value A s w 43.43 (data not shown). It means that, while the 
range of the oscillations increase when k\ goes to zero, at the same time their 
total energy goes to zero. It is understandable: this is because their energy 
goes to zero that their range can increase. On the one hand, it suggests that 
these oscillations poorly contribute to the total energy per unit length of the 
vortex, which should be dominated by the kinetic energy. On the other hand, 
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if fully attributed to the kinetic energy, the above decrease of the oscillations 
energy will be interpreted as an increase of the hollow core radius. 

This paradox was underlined in [8]. Indeed, Equation 33, applied to liquid 
helium at the freezing pressure, gives a deformation energy of 0.0176K/A, 
i.e. within the error bar of the vortex energy (1.90(3)K/A [8]). A precise 
determination of this deformation energy would give an estimate of the roton 
gap, thus completing all the informations on the roton dispersion curve. 

8 Discussion and Conclusion 

We have shown that the density profile of the superfluid close to a singularity 
is expected to be universal when the roton gap Tifl is small. We have derived 
an explicit form of this behavior (Eq. 20) valid as well for real supcrfluids 
as for Gross-Pitaevskii models. We have emphasized that only the additional 
phase-shift </> c and the overall amplitude a depend on the precise shape of 
the model. The linear prediction compares as well with numerical simulation 
of the ground state of the Pomeau-Rica and Bcrloff-Roberts models, as with 
first principle calculations in presence of a vortex. However, while in the first 
case, the order parameter is concerned in a Gross-Pitaevskii approach, in the 
second case, the oscillations directly concern the fluid density. 

This agreement allows to use the oscillations in a reciprocal way, and to 
deduce some characteristics of the roton minimum from their shape. This is 
particularly useful with the first principle calculations where a direct "mea- 
surement" of the dynamic properties of the obtained ground state is very 
difficult [15]. The roton wave number, and the width of the roton minimum 
can be deduced from the wavelength and the decrement of density oscillations 
close to a singularity. Even the roton gap could be deduced from a precise 
determination of the energy of these oscillations. 

As for the vortex, the energy of the corresponding density oscillations is 
negligible compared to the kinetic energy associated with the fluid circulation 
around it. As a consequence, the vortex dynamics is unaffected by these 
density oscillations as far as the curvature of the vortex is smaller than k\ . In 
particular, small amplitude Kelvin waves should be unaffected [16]. However, 
large amplitude Kelvin waves, with wavelength smaller than l/k\ will be 
affected. This point goes out of the scope of this paper. 

While we focused on 4 He close to the freezing pressure, our conclusions 
can be extended to other examples. Metastable overpressurized Helium has 
already been obtained within nanopores [17], or transiently, with acoustic 
oscillations [18]. Many attempts aimed at obtaining metastable superfluid 
molecular hydrogen [19], and estimations give good hope for a next future 
[20] . Finally, even cold atomic gases seem able to sustain roton like excitations 
[21]. In all these cases, the dynamic properties of rotons can be precisely 
deduced from static measurements. 
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